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ABSTRACT 

V407 Vul (RXJ1914.4+2456) and HM Cnc (RXJ0806. 3+1527) are X-ray emitting stars 
with X-ray light curves that are 100% modulated on periods of 569 and 321 seconds 
respectively. These periods are thought possibly to represent the orbital periods of close 
pairs of white dwarfs. In this paper we present optical light curves taken with the high- 
speed CCD camera ULTRACAM on the 4.2m WHT in May 2003 and August 2005 
■ and with the VLT in November 2005. The optical and X-ray light curves of HM Cnc 

' have been reported as being in anti-phase, but we find that in fact the X-rays peak 

around 0.2 cycles after the maximum of the optical light, as seen also in V407 Vul. 
Q^. The X-ray/optical phase shifts are well explained under the accreting models of the 

I ' systems if most of the optical modulation comes from the heated faces of the mass 

O ^ donors and if the X-ray emitting spots are positioned in advance of the mass donors, 

' as is expected given the angular momentum of the accreting material. Some optical 

^ I emission may also come from the vicinity of the X-ray spot, and we further show 

that this can explain the non-sinusoidal lightcurves of HM Cnc. On the basis of this 
model we constrain the temperature of the heated face of the companion star finding 
k>( \ a bolometric luminosity > lO'^'^ ergss^^ and a distance, d > 1.1 kpc. We can identify 

. no explanation for the X-ray/optical phase-shifts under the intermediate polar and 

' unipolar inductor models of the systems. The only significant difference between the 

two stars is that V407 Vul is observed to have the spectrum of a G star. The variation 
in position on the sky of a blend of a variable and a constant star can be used as a 
measure of their separation, and is sensitive to values well below the limit set by seeing. 
We apply this "pulsation astrometry" to deduce that the G star is separated from the 
variable by about 0.027" and hence plays no role in the variability of V407 Vul. We 
show that light travel time variations could influence the period change in V407 Vul 
if it forms a triple system with the G star. 

Key words: binaries: close- stars: individual: V407 Vul, HM Cnc - white dwarfs - 
stars: magnetic fields - X-rays: stars - astrometry 



INTRODUCTION 



V407 Vul (RXJ1914.4+2456 iMotch et all 
and HM Cnc ( RXJ0 806.3-H527, 



Israel et al.l 



199€) 



199C 



iBurwitz fc ReinschI [2003) were both discovered in the 
ROSAT all sky survey and have very simil ar X-ray prop- 
erties . They have period s of P = 569 sec (ICropper et al 



1 19981 : iMotch et all Il996l ') and P = 321 sec jlsrael et al 



[1999') respectively. In each star , only one period (and its 
harmonics) has b een observed (jRamsay et al.l bool I2OO2I : 
llsrael et al] I2OO2I ) at all wavelengths. Taken together, the 
observations have lead to a belief that the periods may 
be orbital, making these the shortest period binary stars 
known, and probably composed of pairs of white dwarfs. 
This would make these systems strong emitters of gravita- 
tional waves and possible progenitors or representatives of 
semi-detached AM CVn stars. 
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There are several rival models for these stars, all of them 
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based upon binary systems. The intermediate polar (IP) 
model llMotch et all Il996l : llsrael et al.l 1 19991 : iNorton et all 
I2OO4I ) is the only one in which the pulsation periods are 
not assumed to be orbital. In this model, the pulsations 
are ascribed to the spin of white dwarfs accreting from 
non-degenerate secondary stars; the orbital periods are pre- 
sumed undetectable. The other three models all invoke 
double white dwarf binaries in which the pulsation peri- 
ods are the orbital periods. There is one detac hed model 
(i.e n on-accreting), the unipolar inductor model (|Wu et al.l 
|2002| ). also called the electric star model because it is pow- 
ered by the dissipation of electric currents induced by an 
asynchronism between the spin period of a magnetic white 
dwarf and the orbital period within a detached double 
white dwarf binary. The other two models each employ 
semi-detached, accreting double white dwar fs: one is mag- 
netic , the double degenerate polar mode l (ICropper et al.l 
1 19981 : iRamsav et "al]l2002l : llsrael et al.ll2002l ) while the other 
is non-magnetic, the direct i mpact model j Nelemans et al.l 
I2OOII : iMarsh fc Steeghsll2002l : iRamsav et al.lf2002h . in which, 
due to the compact dimensions of these systems, the mass 
transfer streams crash directly onto the accreting white 
dwarfs. 

It has proved hard to decide which, if any, of the 
models is correct. Compared to typical accreting systems, 
HM Cue has a weak optical line emission, while V407 
Vul has none at all. This favours the unipolar induc- 
tor model which is the only one without accretion. The 
unipolar inductor model, along with the IP model, is also 
favoured by the observed decrease in pulsa t ion periods 
llStrohmaveij '2002. 2003; iHakala et all l2003l : IStrohmaveil 
I2OO3I : Hakala et al.. .2004 ) although recently accreting mod- 
els with long-lasting spin-up phases have be en developed 
ijP'Antona et al.ll2006l : iDelove fc Taamll2006l) . The shapes 
and phases of the X-ray light curves on t he other hand count 
against the unipolar inductor model (|Barros et al.l |2005| ) 
which can only accommodate the high X-ray luminosity of 
V407 Vul with a white dwarf that spins faster th an its orbit 
l|Marsh fc Nelemansll2005l : iDall'Qsso et"aLll2006al lbl) . The ac- 
creting double-degenerate models on the other hand lead to 
high accretion rates and strong heating of the white dwarf, 
particularly in the case of HM Cnc, which is required to be 
at a d istance of 4 to 20kpc, and well out of the Galactic 
plane (|Bildsten et al.ll2006l : [P'Antona et aLll2006l ). At the 
moment therefore, there is no clear winner, or even lead- 
ing contender amongst the models and better observational 
constraints are a priority. 

Previous studies of the systems have focused mainly 
upon the properties of the X-ray light curves with optical 
data used mainly to track the decreasing periods with less 
attention bei ng paid to the sh a pes of the li ght curves. With 
the w ork of iD'Antona et al.l (|2006l ) and iDelove fc TaamI 
l|2006t ) adding uncertainty to the interpretation of the period 
change measurements, the light curves themselves take on 
more significance. In this paper we present high-speed pho- 
tometry of these systems in three simultaneous bands taken 
in the hope of using the optical characteristics to learn more 
about the systems. In section 2 we report our observations 
and data reduction. In section 3 we present our results. In 
section 4 we use our results to try to determine the origin 
of the optical pulses and explore the consequences for the 
accretion geometry in these systems. 



Target 


Date 


UT 


Seeing, clouds 


V4n7 V11I 


21 


May 2003 


05:33 - 


06:25 


1.0 clear 


V4n7 Vul 


22 


May 2003 


03:28 - 


04:24 


1.0 clear 


V4.n7 Vul 

V t:VJ I V Lil 


22 


May 2003 


04:54 - 


06:25 




V4n7 Vul 

V rtLf 1 V Lll 


23 


May 2003 


02:25 - 


04:24 


1 n f'\f^^v 


V4n7 Vul 

V rtU ( V Lll 


24 


May 2003 


02:48 - 


03:41 


1 n GnTTIP 
J-.U, OUllll^ 


V4n7 Vul 


24 


May 2003 


04:50 - 


06:18 


1 11 /^l no v 
J-.u, cicai 


V4n7 Vul 

V ( V 111 


25 


May 2003 


01:45 - 


02:29 


1 9 plojir 


V4n7 Vul 

V rtLf t V Lll 


25 


May 2003 


03:19 - 


04:41 


1 9 f'\l^^V 
1-.^, L.lV-><H 


V4.n7 Viil 

V '-t\j 1 V U-1 


27 


Aug 2005 


21:10 - 


01:02 


1 1 /^l V 

j_ . 1 , Clival 


V407 Vul 


28 


Aug 2005 


21:05 - 


22:38 


0.9, clear 


V407 Vul 


30 


Aug 2005 


20:50 - 


23:55 


0.8, dusty 


V407 Vul 


31 


Aug 2005 


20:49 - 


22:56 


0.7, dusty 


V407 Vul 


01 


Sep 2005 


20:45 - 


22:58 


0.9, dusty 


HM Cnc 


21 


May 2003 


22:11 - 


23:30 


1.2, clear 


HM Cnc 


22 


May 2003 


21:54 - 


22:57 


1.0, clear 


HM Cnc 


23 


May 2003 


21:57 - 


22:54 


1.0, clear 


HM Cnc 


25 


May 2003 


21:55 - 


22:39 


1.3, clear 


HM Cnc 


27 


Nov 2005 


05:03 - 


06:51 


1.3, clear 


HM Cnc 


28 


Nov 2005 


05:10 - 


08:47 


1.0, clear 


HM Cnc 


29 


Nov 2005 


05:35 - 


08:51 


0.8, clear 



Table 1. Observation log. 



2 OBSERVATIONS AND REDUCTION 

We observed with the hi gh-speed CCD camera ULTRACAM 
iDhillon fc MarshllioOlh mounted on the 4.2m WiUiam Her- 
shel telescope (WHT) in La Palma on May 2003 and Au- 
gust 2005, and mounted on the UTS unit (Melipal) of the 
Very Large Telescope (VLT) in Chile in November 2005. For 
V407 Vul we have observations on five consecutive nights 
from the 21"* to 25"" of May, 2003 with a total of approxi- 
mately 3600 frames of 9.7 sec exposure in the i', g' and u' 
filters and another 2000 frames of 15 sec exposure in five 
extra nights from the 27**" of August to the I''* of Septem- 
ber 2005 in r' , g' and u . For HM Cnc we have around 2000 
frames taken in four nights from the 21''' to 25"" of May 
with 10.1 sec exposures in i' , g' and u' and another 18,000 
frames taken in November 2005 in r', g' and u' with expo- 
sures of 1 to 6 sec. The observing conditions are summarised 
in Table [T] All the times were transformed to TDB, and then 
shifted to time as observed at the solar system barycentre 
using the IDL routine barycen and recorded as a modified 
Juhan day MJD(TDB). 

The data were reduced using th e ULTRACAM pipeline. 
We tried "optimal" photometry (|Navloij Il998h . variable 
aperture photometry and fixed aperture photometry to ex- 
tract the light curves. Optimal photometry gave the higher 
signal-to-noise with the only exception the r' band in the 
August 2005 data, for which we used a fixed aperture ra- 
dius. Optimal photometry requires the profiles to be iden- 
tical in shape and can cause difficulties if this is not the 
case and we believe that in this one case this outweighed 
the improvement in stochastic noise. The subsequent data 
analysis was carried out with IDL. V407 Vul is in a crowded 
field so care was taken to prevent the sky annulus from be- 
ing contaminated by other stars. It is trickier to allow for 
the faint stars that can contaminate the target aperture in 
poor seeing. These are a particular problem in the i' filter 
(May 2003 data) where we found the flux could increase by 
as much as 5% in the poorest seeing. Although relatively 
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few of the data were affected by this, we corrected for it 
by fitting and removing the trend of flux versus seeing from 
the i' data. Finally, the g' data from the second half of the 
second run of the 22"^* of May 2003 and the second half of 
the second run of the 24*'' of May 2003 could not be used 
because V407 Vul was unfortunately positioned close to a 
column of poor charge transfer on the g' CCD. 

In the May 2003 observations of V407 Vul we used two 
comparison stars, one (cl) for the i' and g' bands and the 
other (c2) for the u' images (because cl was too faint in 
It'). The position relative to the target and the magnitudes 
of these comparison stars and the one used for HM Cnc 
are given in Table |3l In the August 2005 observations of 
V407 Vul we only used comparison star c2 because cl was 
saturated in r' due to the longer exposure time. This run 
also suffered from Saharan dust that lead to an extra and 
variable extinction of ~ 0.5 magnitudes at the zenith making 
it impossible to derive an absolute calibration to better then 
0.2 magnitudes. Therefore we used the g' and u' magnitudes 
of c2 calculated in May 2003 to calibrate the August 2005 
data. To obtain the r' magnitude of c2 we applied the same 
correction as for g' . 

The measured mean magnitudes of the systems are 
given in Table [2] As far as possible, the magnitude calibra- 
tion was carried out by comparing the target and the com- 
parison at the same airmass as we did not have sufficiently 
long runs to estimate accurate extinction coefficients. The 
uncertainties of the comparison star for HM Cnc are domi- 
nated by the uncertainties in the extinction coefficients for 
the night both in May 2003 (i) and in August 2005 (r, g 
and u) because in this case we did not observe the target 
and the comparison at exactly the same airmass and some 
correction was needed. 



3 RESULTS 
3.1 Ephemerides 

To compare our optical data with the published X-ray data 
we had to fold our data, on the X-ray epheme ris. Unfortu- 
nately none of the ephemerides p ublished so far ( Strohmavor 
I2004l . l2005l : llsrael et al.]l200l[2004 ,: Ramsay et al. 2006 ) give 
the covariance terms of the fitted coefficients which are 
needed for a correct evaluation of the uncertainties. There- 
fore w e ha d to digitise an d fit th e data of [Strohmaye^ ()2004l . 
I2OO5I ) and iRamsav et al.l (|2006l ) so that we could obtain a 
timing solution whose uncertainties we could compare with 
our data. When we did this we realised that there was an in- 
consistency between the e phemerides of V407 Vul published 
bv lstrohmave r ( 2 0o3 ') and lRamsav et a l. (2006). After inves- 
tigation we concluded that IStrohmavej 's (j2004h ephemeris 
is probably in error because the ROSAT times were not 
corrected from UTC to TT. We therefore use our fitted 
iRamsav et "ail 's l|2006l ) ephemeris (Tabic [XlJ for V4q7 Vu l 
which is similar to ephemeris given in .Ramsav et al.l (j2006l ) 
but has a slightly different z>. For HM Cnc we used 
IStrohmaveij 's (|2005t ) ephemeris. Both ephemerides and re- 
spective covariance terms are given in Appendix A where 
we provide full details of our investigations. 



Filter 


Semi-amplitude 


t-to 







(%) 


(s) 


(cycles) 


i' 


3.03±0.06 


-4.9±1.8 


0.9612±0.0032 


g' 


8.47±0.09 


0.0±0.9 


0.9698±0.0016 


u' 


20.50±0.61 


3.9±2.7 


0.9767±0.0047 


r' 


4.39±0.06 


-0.7±1.2 


0.9596±0.0021 


g' 


8.70±0.07 


0.0±0.7 


0.9607±0.0013 


u' 


21.64±0.44 


O.Oil.S 


0.9607±0.0033 



Table 5. Results of single harmonic sinusoid fitting for V407 Vul. 
The first three lines show the results for the May 2003 data and 
the last three lines show the results obtained in August 2005. The 
times mark the position of the maximum phases and are To = 
52782.191666 for May 2003 Tq = 53612.9483393 for November 
2005. The phases are relative to the ephemeris of Table IXT] 

3.2 V407 Vul 

We show our phase-folded light curves of V407 Vul folded 
on the ephemeris of Table IXT1 on the left of Fig. [1] The two 
datasets (May 2003 and August 2005) were rebinned to 100 
phase bins using inverse- variance weighting to maximise the 

signal-to-noise ratio. 

We computed the Lomb-Scargle periodogram (iLombl 
1 19761 : IScargleiri982l ') and confirmed the previously measured 
period of 569 sec. We then tested how close the signal is to 
a perfect sinusoid by fitting a sine wave at the fundamen- 
tal frequency and at the second and third harmonics. The 
third harmonic is consistent with zero. The results for the 
relative amplitude and the phase at maximum of the funda- 
mental (i.e. the "first harmonic") and the second harmonic 
are shown in Table |4] which also shows the corresponding re- 
sults for HM Cnc. We also fitted a sinusoid with frequency 
fixed to the value derived from the ephemeris of Table lAll 
at our observing date to obtain a normalised amplitude of 
variation and the time (or equivalently the phase) of the 
maximum. The normalised amplitudes, the phase and the 
time-shifts relative to the g'-band are presented in Table [5] 
The amplitude increases strongly towards short wavelengths 
but there is no observable phase shift with wavelength. From 
Table [5] we calculate a difference of phase between our two 
runs of 0.0089 ±0.002. This could be taken to be as a signifi- 
cant shift in phase, however the uncertainty only represents 
the measurement error. When we include the uncertainty 
of the ephemeris calculated with Equation IA2I we obtain 
0.0089 ± 0.019, and therefore we conclude that there is no 
significant variation of the phase shift between the optical 
and the X-rays between the two epochs of our observations 
and that the new ephemeris can be used to extrapolate to 
later times. To compare the optical phases with the X-ray 
light curves it is important to notice that the absolute error 
of the phase due to the ephemeris of Table IXT1 is 0.0090 for 
May 2003 and 0.019 for August 2005. 

3.3 Pulsation astrometry of V407 Vul 

A totally unexpected feature of V407 Vul is that its op- 
tical spectrum is dominated by that of a late G/early K 
star which cannot fit within a 10-minute period binary orbit 
l|Steeghs et al.l [20061 ') . Although this immediately suggests 
the IP model in which one expects a main-sequence sec- 
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Target i' r' g' u' 

V407 Vul May 2003 18.95±0.05 - 20.30±0.06 21.56±0.10 

V407 Vul Aug 2005 - 19.3±0.1 20.29±0.06 21.53±0.08 

HM Cnc May 2003 21.5±0.1 - 20.9±0.1 20.5±0.1 

HM Cnc Nov 2005 - 21.21±0.10 20.77±0.11 20.51±0.12 



Table 2. Magnitudes measured for the two targets. 



Comparison 


Aq arcsec 


AS arcsec 


i' 


r' 


9' 


u' 


V407 Vul cl 


+3.1 


-8.4 


14.21+0.01 


saturated 


16.26+0.01 


19.78+0.03 


V407 Vul c2 


+39.5 


-37.0 


15.73+0.01 


16.08+0.1 


16.96+0.01 


18.84+0.03 


HM Cnc 


-16.9 


-16.4 


15.25+0.11 


15.31+0.10 


16.00+0.11 


17.73+0.12 



Table 3. Positions relative to the target and magnitude of the comparison stars used to flux calibrate the data. 



ondary star (|Norton et al. 1 12004'), the star shows no radial 
velocity variations at a level that rules out o rbital periods 
typic al of cataclysmic variable stars (< Idav. Isteeghs et al.l 
I2OO6I ). Alternatives are that this star is a line-of-sight coin- 
cidence (the field is a crowded one), or it could be part of a 
triple system with the variable. To discriminate between the 
latter possibilities we searched for variations in the position 
of V407 Vul on its 569 period. The idea behind this "pulsa- 
tion astrometry" is that although we cannot spatially resolve 
the variable and G star components of V407 Vul directly, we 
can use the pulsations of the variable to try to detect their 
separation because their mean position will move back and 
forth between the variable and the G star as the variable 
brightens and fades. This method is sensitive to separations 
well below the seeing. 

We measured the position of V407 Vul relative to nearby 
stars in the field. We then computed the amplitude of the 
best-fitting sinusoid over a range of frequencies for both the 
X- and y-positions in each of the three filters as shown in 
Fig. [21 We computed false alarm probabilities using Monte 
Carlo simulatio ns (finding values that agree with the theo- 
retical values of Schw arzenberg-Czernvl [l998l ) . In Fig. [2] we 
show the 99.9% detection threshold for a known period (hor- 
izontal solid lines) and also the 99.9% detection threshold 
for an arbitrary unknown period (dashed lines). We choose 
the 99.9% level because it corresponds to a detection limit 
of about "3a". Note that the detection criterion is more 
stringent when we don't know the period because a penalty 
must be paid for search ing multiple independent periods 
IIHorne fc Baliunas|[l986l ). We include this level to show that 
there are no such detections of any other periodicities. In the 
case of V407 Vul we know the period that we are looking 
for so it is the lower threshold represented by the solid lines 
that applies. As mentioned above, the g' data of the May 
2003 run were partially affected by poor charge transfer in a 
column close to V407 Vul. This has more of an effect upon 
position (especially at the levels we measure here) than on 
flux, so for the position measurements we discarded the 50% 
of the g' data where V407 Vul was closest to the column, 
but as a result reduced the sensitivity of the jr'-band data 
in the left-hand side of the figure. 

There are detections of a signal at the 99.9% level in the 
j/-position data in both r' and g' of the August 2005 run. 



0.05 r 
0.04 hi 
S~ 0.03 ^ 




-0.02 P , , , , I , , , , I , , , , I , , , , j 
0.0 0.5 1.0 1.5 2.0 

Orbital Phase 



Figure 3. Phase folded position variation for g' and r' the 
datasets that show a significant signal. The arrow shows the po- 
sition of the maximum of the flux. 

Fig.|3]shows a phase-folded, binned plot of the {/-position for 
these two cases. The time of maximum excursion roughly 
corresponds with the time of maximum light as expected, 
and both datasets are consistent with each other in this 
respect. However, the signal is tiny, with an amplitude of just 
0.005 pixels or 0.0015 arcsec, and so we endevoured to test 
the reliability of this detection as far as we were able. The 
most obvious problem is that V407 Vul is in a crowded field 
and so the position measurements could be affected by other 
stars. There are two stars w ithin 1.5 arcsec of V4 07 Vul that 
can be seen in Figure 2 of lRamsav eT^ (|2002D . To check 
how these stars affected our measurements we first tested 
whether the detection depended upon the FWHM of the 
seeing. We divided the data in two parts, higher and lower 
FWHM. The reduction of data size lead to no detection 
in either case but the significance of the peaks was higher 
in the small FWHM dataset. The reverse would have been 
expected had blending with the two nearby stars been the 
cause. 

We measured the centroids by cross-correlation with 2D 
Gaussians of fixed width. This allows us to assess the effect 
of the Gaussian width upon the measured amplitude. As the 
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Figure 1. Phase-folded light curves of V407 Vul (left) and HM Cnc (right) using Table IAH and IStrohmaved 's Hooi). The flux is 
normalised to unity in each case. The different filters are displaced vertically for clarity. 



Filter 




V407 Vul 






HM Cnc 






02/11 




(t>2 


a2/ai 


01 


02 


i' 


0.079±0.030 


0.961±0.012 


0.458±0.16 


0.207±0.064 


0.122±0.010 


0.335±0.025 


g' 


0.053±0.014 


0.970±0.005 


0.388±0.08 


0.157±0.025 


0.117±0.004 


0.338±0.012 


u' 


0.095±0.036 


0.977±0.012 


0.544±0.05 


0.131±0.050 


0.118±0.008 


0.285±0.031 


r' 


0.024±0.015 


0.960±0.020 


0.361±0.099 


0.202±0.009 


0.1616±0.0014 


0.362±0.0034 


g' 


0.039±0.012 


0.961±0.013 


0.444±0.14 


0.188±0.005 


0.1659±0.0008 


0.356±0.0022 


u' 


0.012±0.022 


0.961±0.033 


0.340±0.27 


0.205±0.014 


0.1700±0.0022 


0.342±0.0055 



Table 4. First and second harmonic decomposition of the optical light curves for V407 Vul and HM Cnc. ai and 02 are the sem i- 
ampli tudes of the first and second harmonics respectively and (j>i and (j>2 their phases of maximum light on Table ED s and IStrohmaveit s 
1120051 ) ephemeris. In the case of HM Cnc, the measurements at the top and bottom come from the WHT and VLT respectively, hence 
the marked difference in the uncertainties. 



X 



FWHM of the Gaussian increases, we expect to see a more 
pronounced impact of the nearby stars. Therefore if it is the 
nearby stars rather than the G star that are responsible for 
the variation, we expect an increase of measured amplitude 
with Gaussian width. In fact we see the reverse as Fig. |4] 
shows, at least in the y-positions for which we have detec- 
tions. The X positions do show a distinct upturn at large 
FWHM o wing to the much brig hter star 5" East of V407 Vul 
(star B of iRamsav et al.ll2002i ) which was positioned to the 
left of V407 Vul in our data. 

As a final check we carried out simulations of our po- 
sition measurements using parameters matching the stars 
that we could see nearby V407 Vul, including the two very 
close ones mentioned above. This leads to the dotted line in 
Fig. |4l In viewing this figure it must be recognised that the 
data are not independent and so to some extent the trends 
with FWHM can just reflect noise; the dashed lines in the 
figure show two simulations of the effect that noise can have 
upon the simulated amplitudes. These show that for the y- 
positions the measured amplitudes are indeed significantly 
larger than the simulated values, and provide further confi- 
dence in the reality of the detection. 

We conclude, albeit tentatively, that we have detected 
a change in the spatial position of V407 Vul that is corre- 
lated with its pulsations and that the change in position is 



0.005 




234567 234567 
Gaussian FWHM Gaussian FWHM 

Figure 4. The dependence of the amplitude of the variation of 
the position upon the width of the Gaussian used to calculate 
the position. The data are plotted as a solid line, the simulations 
accounting for the known stars are plotted as dotted lines (no 
noise) and dashed lines (with noise). The simulations with noise 
have similar errors to the data but these errors were not plotted 
for clarity. The plate scale is 0.3" per pixel. 



because the G star we see in its spectrum is not exactly coin- 
cident with the variable. We predict that the G star should 
be below V407 Vul in our field which roughly corresponds 
to south of V407 Vul. We obtain amplitudes of the posi- 
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Figure 2. Each panel shows the amplitude spectra of the x (left) and y (right) variation of the position of V407 Vul for the three filters 
i' or r', g' , and u' from top to bottom. The left panel shows the May 2003 data while the right shows the August 2005 data. The vertical 
dashed lines show the position of the 569 sec period. The solid horizontal lines show the 99.9 percent significance level for a known period 
and the dashed horizontal lines show the same level for an arbitrary period. 



tion variation which we denote by p of 0.00512 ± 0.0012 and 
0.00514 ± 0.0010 pixels for g' and r' respectively. 

The value of p is related to the separation on the sky d, 
the fractional amplitude of the flux variation a as listed in 
Table [5] and the fractional contribution of the G star to the 
flux at minimum light /, through the following relation; 



P = 



(1) 



Using the measured values for a and p we calculate f^id — 
0.0366" ± 0.0073 and fg,d = 0.0192" ± 0.0046. This gives a 
value of fgi/fr' = 0.52 ± 0.16. T his is consisten t with the 
spectra of the G star from which ISteeghs et al.l (|2006l ) es- 
timate that fr' > 0.85 and fgi > 0.6. These numbers also 
match the amplitude of the flux variation whose significant 
drop from u' to g' to r' (Table [S} can be explained by di- 
lution of an underlying variable with a constant amplitude 
with wavelength, as for HM Cnc. If we assume fgi — 0.7 we 
obtain d ~ 0.027"; this compa res with the upper limit of 
0.1" set bv lSteeghs et all (|2006|). The distance to the G star 
of 1 kpc estimated by Steeghs et al.l l|2006f ) , leads to a mini- 
mum separation of ~ 30 AU, equivalent to a period of 120 yr, 
and means that the G star cannot be the direct cause of the 
optical and X-ray pulsations. Nevertheless it may well be as- 
sociated with the system in the form of a hierarchical triple, 
a point we return to after we have presented the lightcurves 
of HM Cnc. We finish by noting that our failure to detect 
anything in the u' band is to be expected. Assuming typical 
colours for the G star and a hot spectrum for the variable, 
we expect that if f'g = 0.7, then /4 — 0.3. The effect of this 
reduction in /, which is to make any movement more diffi- 
cult to detect is in large part offset by a factor 2.1 increase 
in a/(l + a), but then we are faced with a factor 8 worse 



sensitivity in the u' band, and the result is that there is no 
detection in the u' band data. 



3.4 HM Cnc 

We present the phas e-folded light curves of HM Cnc us- 
ing jstrohm^eJslloO^ ephemeris in the right-hand panel of 

Fig.m 

We computed the Lomb-Scargle periodogram to confirm 
its 321 sec period and we noticed that the relative strength 
of the second harmonic is high er than V407 V ul's, as has 
already been pointed out bv Jsrael et al.l (|2002l ). This is in- 
deed clear from the non-sinusoidal shape of the light curves 
in Fig. [T] The results of the relative strength of the first and 
second harmonic and their phases are shown in Table [3] The 
second harmonic is approximately 15% of the fundamental 
and its maximum occurs 0.2 of a cycle after the maximum of 
the fundamental. This results in an asymmetry in the light 
curve whose rise time is longer than its decline; we discuss 
its origin in section 3] 

We applied the same method as for V407 Vul to ob- 
tain the normalised amplitude of variation and the time and 
phase of the maximum. These results are presented in Ta- 
ble [S] In this case the amplitude of the variation decreases 
slightly for shorter wavelengths which reinforces the picture 
that in V407 Vul the change with wavelength is due to di- 
lution at long wavelengths by light from the G star. For 
HM Cnc the normalised amplitudes of variation are smaller 
than the u' band for V407 Vul (which is the least contami- 
nated by the constant star). This could be easily explained 
by the inclination of the plane of the orbit and/or differences 
in temperatures of the stars. 

The higher signal-to-noise ratio of the VLT data from 
November 2005 reveals a trend with waveband in the phase 
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Filter 


amplitude 

(%) 


i-to 
(s) 


(cycles) 


i' 


14.77±0.95 


1.7±3.2 


0.121±0.010 


g' 


13.48±0.34 


O.Oil.3 


0.116±0.004 


u' 


13.08±0.66 


0.5±2.6 


0.118±0.008 


r' 


13.54±0.12 


-1.30±0.45 


0.1615±0.0014 


g' 


12.74±0.07 


0.00±0.27 


0.1656±0.0009 


u' 


11.90±0.17 


1.32±0.72 


0.1697±0.0023 



Table 6. Results of single harmonic sinusoid fits to the HM Cnc 
data. The times mark the phases of maximum light and are ref- 
erenced to two times: Tq = 52782.895768 for the May 2003 data 
and To = 53702.3368167 for the November 2005 data. 



of the fundamental which is progressively delayed towards 
short wavelengths. To test whether the trend is significant, 
we carried out an F-ratio test comparing two models, one 
of a constant phase in the three bands versus one of a lin- 
ear trend of phase with wavelength, using the central wave- 
lengths of each band: 3543, 4770 and 6222 A. The F-ratio 
is the ratio between the x'^/ (number of degrees of freedom) 
of one fit divided by the same quantity for the other fit. 
We only had three points so the constant model has two 
degrees of freedom while the straight-line fit has just one. 
The values of the are 10.82 and 0.068 for the constant 
and straight-line respectively, giving an _F-ratio of 79.5. This 
is significant at the 90% level but not at 95%, so, although 
suggestive, there is no significant shift. 

Table[6]shows that there is a phase difference of 0.050± 
0.004 between our two runs (May 2003 to November 2005) 
where this is the measurement error only. As with V407 Vul, 
we also have to add the uncertainty of the ephemeris (see the 
appendix). The error of the difference of phases due to the 
uncertainty of the ephemeris calculated using Equation IA2I 
and the correlation coefficients given in Table 1X21 is 0.013. So 
there is a phase difference between the two runs of 0.050 ± 
0.014. Therefore there is marginally significant variation in 
phase which might mean that there is a variation of the 
phase shift between the optical and the X-rays or, more 
likely, that the spin up rate is varying. The uncertainty in 
the absolute phase calculated using Eguation lAll is 0.005 for 
the May 2003 data and 0.01 for the November 2005 data. 
These are useful to compare the optical phases with the X- 
ray phases, and as we shall see next there is a significant 
phase shift between the two. 



3.5 The Optical/X-Ray phase shift of HM Cnc 

The relative phases of the optical and X-ray light curves 
arc an important constraint upon models. [Israel et al.l l|2003l . 
[2004 ) found that optical and X-ray light curves of HM Cnc 
were in anti-phase as might be expected for an X-ray emis- 
sion region facing the secondary star, contrary, for example, 
to expectations based upon the direct impact model. 

In the right-hand panel of Fig. [ 5] we present t he X - 
ray and optical light curves folded on IStrohmaved' s (120051) 
ephem eris. Our phase shift differs from Israel et ahl 's ( 200^: ; 
[2OOJ) studies by around 0.2 cycles. To test if this was a gen- 
uine change in the system, we reduced some of the archival 
i?-band VLT data from the 12*'' of December 2002 used 



2.5 
2.0 
1.5 
1.0 
0.5 

X 

.30.0 



u 

g 

f?fsi~i, ^iStWs'yj:^'* ^l*i«i«i«-'--f»'m- StM-^stf^tF^ik i^f^ftrinfr 

i 


V407 Vul u 

. _______ g ; 

r 







I 200 400 600 800 



100 200 300 400 500 600 700 
RXJ0806 I 



1000 



-1 

100 200 300 200 400 600 800 
Frame number Frame number 

Figure 6. V407 Vul (top panel) and HM Cnc (bottom panel) 
light curves after removal of the sinusoidal pulsations. For each 
object the data from the first observing period is in the left and 
the second in the right. There are no significant variations of 
either source. We inserted gaps between different nights. 



bv llsrael et al.l (|2004l ). We reduced these data with the UL- 
TRACAM pipeline and applied the same methods and time 
conversions as for the WHT/ULTRACAM data. The resuhs 
are also shown in the right-hand panel of Fig. [S] and agree 
nicely with our ULTRACAM data. Clearly the system phase 
is stable and there are no problems with the times of either 
data set. The difference must be due to the data reduction. 
We confirmed our timing results with three different data re- 
duction packages so we believe that our relativ e phase is cor - 
rect a nd sus pect that there is a problem with 
(2003; 20041 ) values. We were able to confirm 
(20031 : [2OOJ) X-ray phase so assume that there is a problem 
only with the optical timings. The shift of 0.2 cycles is about 
1 minute, which is suggestively close to the 64 seconds off- 
set between UTC and TDB. Dr Israel was kind enough to 
confirm that such an error was possible. 



Israel et ahl's 



Israel et ahl 's 



3.6 Flickering 

The random stochastic variations known as "flickering" are 
one of the hallmarks of accreting systems. We therefore 
looked for any signs of flickering in our data. A plot of the 
light curves after removing the sinusoidal variation is shown 
in Fig. 

The light curves are very constant except for long 
timescale variability of HM Cnc during the November VLT 
run. The observations of HM Cnc started at high airmass, so 
some of variations seen could be a consequence of extinction, 
except that it' does not look much more variable than g' or 
r'. Therefore we believe that this may be true variability of 
the source and not an artefact. 

In both systems the magnitudes measured in the two 
observing runs agree well within the errors. In the case of 
V407 Vul we also searched for any flux variation on longer 
time scales. We had data from the "auxiliary port" of the 
WHT taken on the lO**" of April 2003 and also Liverpool 
Telescope data taken on the 5"^ of September 2004. The 
different data sets were all within 10% of each other. 
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Figure 5. X-Ra y /optical phase folded light curve of V407 Vul (left) and HM Cnc (right) using the ephemeris of Table lAll and 
llStrohmavei]l2005l) respectively. For the X-ray light curves of V407 Vul we re-analysed Chandra observations, with the solid line showing 
data from 19*'' February 2003 and the dashed line from 24"^ November 2003. For HM Cnc, the X-rays were a dapted from | Strohmaverl 
1I2OO5I '). We overplot the optical g' band showing our results. For HM Cnc we also show the VLT/FORS data lllsrael et al.ll2002l) taken 
in the R filter (top). 



We estimated the variability of these systems by cal- 
culating the RMS of the lightcurves after removing the si- 
nusoidal variations. We filtered the short-term variations to 
minimise the photon noise. We also filtered longer term vari- 
ations so we could compare our short runs with a longer run 
on the cataclysmic variable SS Cyg. For the August 2005 
V407 Vul run we obtain an RMS variability of 0.7% in r' 
and 0.8% in g' . For the November run on HM Cnc we ob- 
tain an RMS variability of 1.6% in r' and 1.0% in g' . We use 
these runs as they have the highest signal-to-noise, neverthe- 
less the variability still contains a significant component of 
photon noise. We do not quote the variability in u' because 
it is completely dominated by photon noise. For comparison, 
applying the same filtering of the data to data on the well- 
known CV SS Cyg, we obtain an RMS variability of 3.0% in 
r' and 5.0% in g' . As mentioned above, the fraction of the 
G star in the g' band of V407 Vul is of order 70%. Therefore 
its intrinsic variability is of order 2.5%, assuming the varia- 
tions are not dominated by photon noise. This is a factor of 
two less than SS Cyg and a factor of fou r if one accounts fo r 
SS Cyg's dilution by its secondary star (|North et al.ll2003 ). 
Other cataclysmic variable stars we looked at are similar 
to SS Cyg, so we conclude that the measured variability in 
V407 Vul and HM Cnc is much less than in normal cata- 
clysmic variable stars. 

The lack of obvious fiickering is a point against accret- 
ing models, although not a conclusive one as there are wide 



variations in the amount of fiickering shown by definitively 
accreting binaries, and one cannot be certain that it should 
have been detected. It does however suggest that most of 
the optical fight does not come directly from the accreting 
region. 



4 DISCUSSION 

4.1 The X-ray versus optical phases 

The correction to the relative X-ray versus optical phase 
of HM Cnc that we have identified makes it very similar in 
this respect to V407 Vul: in each system the X-ray flux peaks 
~ 0.2 cycles after the optical flux. This can be added to the 
shapes of the hght curves as evidence that these two stars 
are indeed related systems, as is evident from Fig. [5] We 
will now investigate what this result implies for the different 
models. 

In the majority of models, the X-rays come from a spot 
on the primary star which moves in and out of view as it 
rotates. The exception is the IP model where the modula- 
tion is the result of the accretion stream flipping from one 
pole to the other although it seems unlikely that such a 
process can really switch off the X-rays as completely as 
observed. Less attention has been paid to the optical pulsa- 
tions. Within double degenerate models, these seem likely to 
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originate from the heated face of the secondary star which 
would naturally explain their near-sinusoidal shape and, per- 
haps, the absence of flickering. Such heating may be a result 
of the X-ray emission from the primary star, or the primary 
star could simply be hot as a result of compressional heating 
(|Bildsten et al.H2006D . 

Assuming that we are correct about the main site of 
optical emission. Fig. [7] shows the geometrical arrangement 
that explains the relative phases of the optical and the X- 
ray light curves. In this model the 0.2 cycles delay of the 
X-ray peak relative to the optical peak implies that the X- 
ray emission spot is rotated ~ 0.3 cycles from directly facing 
the secondary star, in the direction of the orbit. 

This is not the whole story however, because in HM Cnc 
at least the optical light-curve is somewhat saw-toothed in 
shape. As our Fourier decomposition shows, this is caused 
by a significant second harmonic that happens to peak at 
the same phase as the X-rays, as seen in the difference 
(f)2 ~ (pi ~ 0.2 in Table |4l The natural explanation for this 
is that the X-ray emission spot is also the site of some op- 
tical light. If this is a localised region so that the shape of 
its light curve can be approximated by a truncated sinusoid 
{f{e) = cose for -7r/2 < 6 < 71/2, f{e) = otherwise), 
then it can produce a second harmonic. It will also con- 
tribute some first harmonic as well, which means that the 
first harmonic emission that we see is the combination of 
contributions from the heated face of the secondary star and 
the spot on the primary star. This retards the optical phase 
so that the observed 0.2 cycle shift is an under-estimate of 
the true shift between the emission from the heated face and 
the X-ray emission. 

If the optical emission truly can be approximated by the 
truncated sinusoid, then for HM Cnc we find that we can fit 
the phases and harmonic amplitudes if the X-ray spot leads 
the heated face of the secondary star by ~ 0.26 cycles (i.e. a 
little more than 90°) and the optical emission from the spot 
on the primary star peaks at ~ 75% of the amplitude of the 
emission from the heated face. While this is probably rather 
simplistic, it demonstrates that the simple model illustrated 
in Fig. [7] is capable of explaining some secondary details of 
the data. With this decomposition of the optical light, the 
X-ray emission site is ~ 90° ahead of the secondary star 
and it is then not clear if the X-rays can directly heat the 
secondary star or not. V407 Vul has a much weaker second 
harmonic, and so in this case the spot is presumably the full 
0.3 cycles or ~ 110° ahead of the secondary star and cannot 
see it directly. 

The X-ray/optical phase-shifts in V407 Vul and 
HM Cnc are very naturally explained by both the direct 
impact and double-degenerate polar models. The accretion 
spot in normal polars is obser ved to lead the secondary star 
by of order 0.1 to 0.3 cycles l|Croppeij|l988l ) and a similar 
shift is expected in the direct impact mode l, depending upon 
the system parameters (|Marsh et al.ll2003 ). As Fig [7] shows, 
this is exactly what is required to match the observations. 
On the other hand, as far as we can see, there is no natu- 
ral explanation for the phase-shift in the unipolar inductor 
model for which one would expect anti-phasing, unless there 
is some as yet undiscovered mechanism for displacing the 
magnetic footprint of the secondary star in advance of its 
orbit. This is difficult given that the orientation of the pri- 
mary star relative to the secondary star changes relatively 



rapidly in the unipolar inductor model and so a fixed orien- 
tation is hard to contrive. The X- ray /o ptical phase-shift is 
also a difficulty for [Norton et all 's l|2004l ) IP model for which 
they also predict anti-phasing with the optical pulses ap- 
pearing as the accretion stream switches to the hidden pole, 
the X-rays going to zero at this point. IPs are sufficiently 
complex that an offset as observed would not perhaps be 
that surprising, but in any cas e there are other mor e serious 
difficulties with the IP model l|Cropper et al.ll2003 '). 

The direct impact model can be used to predict the 
phase shape and thus, if it is true, constrain the masses of 
the binary components. We define the impact angle as the 
angle between the X-ray emission site and the secondary 
star. We calculate the impact angle for over a grid of Mi 
and AI2. In Figure |5] we show contours of same impact angle 
where the shaded areas represent the probable regions where 
the systems lie. The uncertainties in the angle are higher 
for HM Cnc than for V407 Vul because of the existence 
of the second harmonic component. In the same figure we 
also plotted the dynamic stability limit (dashed line). We 
conclude that for V407 Vul 0.4 < Mi < 0.55 Mq and 0.08 < 
M2 < 0.4 Mq and for HM Cnc 0.6 < Mi < 0.9 Mq and 
0.12 < M2 < 0.45 Mq. 

4.2 A limit on the bolometric luminosity of 
HM Cnc 

Pursuing the idea of the heated face further leads to a lower 
limit upon the bolometric luminosity of HM Cnc, assum- 
ing that the double degenerate models are correct (unfor- 
tunately the G star in V407 Vul's spectrum precludes the 
same calculation). The idea is to derive a lower limit on 
the temperature of the heated face from the spectrum of 
the pulsations, which since it is a measure of the fiux from 
the primary star at a distance equal to the orbital separa- 
tion, which is approximately known, gives us a luminosity. A 
slight complication is that we do not know for sure whether 
the X-rays or the primary star's photosphere is responsible 
for the irradiation. This ultimately leads us to two different 
possible limits. We begin by obtaining the weaker of the two 
which applies in the case of X-ray heating, and then consider 
the revised limit necessary if the photosphere is responsible 
for the heating. 

4-2.1 Temperature of the heated face 

We first derive a lower limit on the temperature of the 
spectrum using our magnitudes and thos e reported by 
iRamsav etlii] (|2002l ). [Israel et al.l l|2002l ) and lReinsch et all 
(|2004l '). as shown in the left-hand panel of Fig. |9l Two black- 
body spectra are shown, each scaled to give the minimum x^- 
One, with temperature 32,400 K (solid line) is the global best 
fit, while the other (dotted), with temperature 18,500 K, has 
the minimum temperature that gives a within the 99% 
confidence threshold. We take this to be the minimum pos- 
sible temperature of HM Cnc given its optical and infra-red 
fluxes. We assume further that this reflects the temperature 
of the primary star, Ti , since if it is the secondary star, the 
primary star would have to be extremely hot to produce 
a significant reflection effect; there is no equivalent upper 
limit as the optical and IR fluxes do not constrain the max- 
imum temperature at all. Armed with the lower limit of 
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Figure 7. In the left-hand panel we show a schematic picture of the primary star (left) with an X-ray emitting spot which has a peak 
of X-ray emission in the direction indicated by the dashed arrow. The relative sizes of the two stars are drawn to match masses of 
Ml = 0.53 Mq and M2 = 0.12 Mq (see text). The secondary star, which orbits counter-clockwise, has a heated face (shaded) whose peak 
emission is in the direction of the solid arrow. The figure is arranged to give the optical (solid) and X-ray (dashed) light curves shown 
on the right, which have the same relative phasing as both V407 Vul and HM Cnc. 
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Figure 8. We show the system mass constraints for V407 Vul (left-hand panel) and HM Cnc (right-hand panel) if we assume the direct 
impact model. The upper dashed line shows the dynamic stability limit. We show contours of equal impact angle. Note that the maximum 
impact angle is approximately 130° which corresponds to the transition between the disc and direct-impact accretion. 



Ti > 18,500 K, we can then use the pulsation amplitudes 
to place a similar lower limit upon the temperature of the 
irradiated face of the secondary star Tirr using the ratio of 
black-body spectra, as shown on the right of Fig. [9] Again 
taking the 99% confidence limit, we find that the tempera- 
ture of the irradiated face must be at least Tirr > 14,800 K; 
this limit rises to 21,700 K if we use the best-fit value for T\ 
of 32,400 K. The lower limit on the temperature of the irra- 
diated face leads directly to a lower limit on the bolometric 
luminosity of the primary star since assuming that the irra- 
diation dominates the intrinsic luminosity of the secondary 
star we have 

iboi = 4-7ra^crT;rr, (2) 

where a is the separation and a is the Stefan-Boltzmann con- 
stant. The strictest lower limit comes from taking the small- 
est separation, which corresponds to the smallest masses 
for the two component stars. We used M\ = 0.53 Mq and 
M2 = 0.12 Mq which ensure that the secondary star can 



fit within its Roch e lobe and that mass transfer is stable 
(|Marsh et al.ll2004h : a smaller value for M\ could be used 
if the system is detached, but would be largely compen- 
sated for by the need for a higher value of M2 to avoid mass 
transfer. Our masses and the period of 321 seconds imply 
a separation of a = 0.089 R0 (these values were used to 
scale Fig. (Tjl. Scaling from the Sun we therefore find that 
Lboi > 0.34 Lq = 1.3 X 10^^ergss"\ 

This is already a significant lower limit as it is some- 
what higher than, but consistent with, the X-ray luminosity 
of Ljf_~_5 X 10^ ergs s~^ at 500 pc distance estimated by 
Isra el et all (|2003l ). However, we have pushed the tempera- 
ture to marginally acceptable values. For instance, the best- 
fit temperature Ti = 32,400 K which leads to Tirr > 21,700 K 
raises the luminosity l imit by a factor of 4.4, which hints at 
a larger distance than llsrael et al.] l|2003l ) assumed. 

Comp aring with white dw arfs of similar temperature 
and mass (|Bragaglia et al.ll 19951 ). the absolute magnitude of 
the primary star is bounded by My < 10.7. Given V — 21.1 
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Figure 9. In the left-hand panel we show the mean fluxes of HM Cnc compared to two black-body spectra. One is the best fit (32,400 K, 
solid), the other has the lowest temperature consistent with the data (Ti = 18,500 K, dotted). The right-hand panel shows the pulsation 
amplitudes that we measure compared to the ratio of two black-bodies (BxiTiy-^) / Bx{Ti)) to place a lower limit upon the temperature 
of the pulsed light T[^^ > 14,800 K; the dashed line shows the limiting case. 



l|lsrael et al.ll2002l ). and assuming that reddening is negligi- 
ble, this suggests that d > 1.1 kpc, although this limit can 
be lowered if we adopt a higher mass for the primary star. 

These limits apply if it is the X-rays that drive the 
heating, but it may well be that it is the photosphere of 
the primary star itself that is important. The spot position 
~ 90° ahead of the secondary means that the spot may not 
be able to see the secondary star at all. This leads to a simple 
but important modification as we now show. 

4-2.2 Heating driven by the primary star's photosphere 

If the primary star's photosphere drives the heating then 
this sets a relation between Ti and Tirr 



(3) 



where T2 is the temperature of the unhealed photosphere of 
the secondary star, Ri is the radius of the primary star, a is 
the orbital separation and we have assumed that all incident 
flux is absorbed and for simplicity we do not try to account 
for the range of incident angles over the heated face, but 
just consider the sub-stellar point. As we said earlier, if T2 
is significant, it is hard to get much of a refiection effect, so 
we take it to be negligible and therefore 

r.^(-) T... (4) 

The masses adopted above give the smallest ratio of a/Ri 
leading to Ti = 2.58Tirr. This equation was used to boot- 
strap from the lower limit on Tirr to a new lower limit of 
Ti, which was then used to place a new lower limit on Tirr 
using the procedure of the previous section. We obtained up- 
dated limits as follows: Trr > 34,800 K and Ti > 90,000 K. 
The heated face temperature rises by a factor of 2.4, and 
so the lower limit on the bolometric luminosity rises by a 
factor of 33 to L > 10 Lq = 4.0 x 10^* ergs s"\ Again 
comp aring with white dwa rfs of similar temperature and 
mass l|Bragagha et al.lll995l ). the absolute magnitude of the 
primary star is bounded by Mv < 8.0. Given V — 21.1 



llsrael et all l|2002l ). and assuming that reddening is negligi- 
ble, we must have d > 4.2 kpc. This would place HM Cnc 
more than 2.5 kpc out of the plane, and it would possibly be 
a halo object. We note that a halo-like transverse velocity 
of 200kms~^ and our distance limit imply a proper motion 
< 0.01 arcsecon ds/yr, below the lim it of 0.02 arcseconds per 
years placed bv llsrael et all l|2002h . 

Our distance limits do not discriminate between accre- 
tion m odels which work best for large distances, in excess of 
4 kpc teildste n et"al] I2OO6I : Ib'Antona et"aLll2006h and the 
unipolar inductor mod el which works well for d < 1 kpc 
iDairOsso et al.l l2006al lbh. However, they do suggest that 
UV observations may have a value in tightening the lower 
limits upon temperatures and hence the distance. 



4.3 Direct impact or polar? 

We have lumped the accreting double-degenerate models, di- 
rect impact and polar together as "accretion models" , as we 
think they provide equally good explanations for our data. 
Fo r V407 Vul the double-degenerate polar model suggested 
bv lCropper et aL IIT998il w as disc arded when no polarisation 
was found llRamsav et al.ll2002l ). In the case of HM Cnc, 
iReinsch et al. I l|2004l ) have claimed a detection of circular 
polarisation but at a low level given the faintness of the tar- 
get (0.5%) that needs confirmation. However, we think that 
the polar model may have been written off prematurely as 
there are some very high-fi eld polars which sh ow very lit- 
tle polarisatio n (AR UMa, [S chmidt e t al.lll996l : V884 Her, 
Szkodv ct al. 1995I . ISchmidt et al. 20oi)and strong soft X- 
ray components, very much like V407 Vul and HM Cnc. 
It has been suggested that this is because the shocks are 
buried in these systems, due to the high accretion rates, 
rather as iMarsh fc Steegt3 (|2002l ) suggested for the direct 
impact model. Polars show stronger optical line emission 
than either V407 Vul or HM Cnc, but this is not a strong 
argument against the polar model since the systems, if they 
are double degenerates, would be helium-rich and very com- 
pact, and so different from normal CVs. 
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4.4 Period changes in V407 Vul and the G star 

We have shown that the G star does not play a direct role 
in the variability of V407 Vul but it could be gravitationally 
bound to the variable, in which case it may cause an appar- 
ent period change through variable light travel time effects. 
How significant could this be? Assuming that the G star 
has mass M, then the maximum acceleration of the binary 
along the line of sight is ~ GcM/a^ where a is the separa- 
tion of the binary and the G star. The subscript c in the 
gravitational constant is to avoid confusion with the G star. 
This leads to a quadratic term in the usual To + PE + CE^ 
ephemeris equal to 



C 



2ca2 



(5) 



where c is the speed of light and P the orbital period. This 
leads to an apparent rate of period change given by 



P = 



GMP 



(6) 



Taking a to be comparable to the projected separation at 
Ikpc of 27 AU, and M = IMq gives |P|max ~ 1.6 x 
10~^^ s/s. This is about 5 times larger than the observed 
value l|Strohmaveij|2005l ) and thus we conclude that the G 
star has the potential to have a significant effect upon the 
rate of period change measured in this system. This adds an 
extra uncertainty that may allow both the unipola r inductor 
l|Marsh fc Neleinansll2005l: iDaH'Osso et al.ll2006al lbh and ac- 
creting models ( D'Antona et al.ll2006l ) to match this system 
more easily than has been the case to date. Continued ob- 
servations in order to place limits upon or detect a relative 
proper motion between the variable and the G star would 
be of interest for testing the triple star model. We estimate 
the orbital velocity of the G star to be ~ 3kms~^, which is 
perhaps detectable given a long enough period of time. 



5 CONCLUSION 

We have presented optical photometry of V407 Vul and 
HM Cnc in i' , r' , g' and u' bands taken with the high-speed 
CCD camera ULTRACAM. For V407 Vul we have a hint 
of detection of a third component in the system at 0.027" 
from the variable. We believe this to be the G star that is 
seen in the spectrum of V407 Vul, which therefore cannot 
be the secondary star of the variable. We cannot distinguish 
whether it is a line-of-sight coincidence or a triple system. 

For HM Cnc we find a new phasing of the X-ray and op- 
tical data which renders it indistinguishable from V407 Vul 
with the optical pulses 0.2 cycles ahead of the X-ray pulses. 
The offsets are naturally produced by double-degenerate ac- 
creting models of the systems, both polar and direct impact, 
but seem hard to reconcile with the unipolar inductor and in- 
termediate polar models. The optical light curves of HM Cnc 
are non-sinusoidal and a Fourier decomposition shows that 
there is likely a contribution to the optical light from the 
same site as produces the X-rays. 

On the assumption that the optical pulses of HM Cnc 
are the result of irradiation of the secondary star within a 
double degenerate binary, and using the relative constancy 
of the fractional pulsation amplitude with wavelength, we 
place a lower limit on the distance to the system of > 1.1 kpc. 



If it is the photosphere of the accretor rather than the X-ray 
site that is responsible for the heating, then this limit rises 
to d > 4.2 kpc. Space ultraviolet observations are the best 
hope for strengthening these constraints. 

Finally we remark that both the polar and direct im- 
pact models provide equally good explanations of our obser- 
vations and that there are high magnetic field polars that 
show similar properties to V407 Vul and HM Cnc i.e. very 
soft X-ray spectra and low polarisation. 
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to (TDB) 

;>(Hzs-i) 
r{<l>o, f) 

r(i/, u) 



49257.533373137 
0.003(30) 

0.00175624626(39) 
9.9(1.9) X 10-1** 
-0.92074289 
0.86174740 
-0.98817908 



Table Al. Epheme ris of V407 Vul derived from the data of 
iRamsav et al.l l|2006l).The uncertainties of the parameters are 
given within parentheses. We also give the correlation coefficients 
for the fitted parameters. 



APPENDIX A: EPHEMERIDES 

AO.l V407 Vul's ephemeris 

As we mentioned in Section 13.11 in order to calculate the 
uncertainty in the published ephemerides we need the co- 
variance terms of the fitted coefficients that are not given 
in any published work. Therefore we had to digitise and fit 
the X-ray data in order to obtain a timing solution whose 
uncertainties we could compare with our data. To digitise 
the data we applied the Linux utilities pstoedit and xfig to 
the PostScript figures from the published papers to obtain 
the coordinates of the points and their error bars. Such a 
process can at best match the original data, and can po- 
tentially degrade it, but in this case the precision of the 
PostScript data is good enough that it has no measur- 
able effect; we were able to confirm our numbers directly 
in one case after Dr Ramsay kindly sent us his data. For 
V407 Vul's ephemeris w e digiti sed the bottom panel of Fig- 
ure 6 from IStrohmaverl (|2004|) and Figure 1 from the re- 
cently published ephemeris of Ramsav et al.l l|2006l ) that ex- 
tends the ephemeris to April 2004. These figures show the 
residuals in phase relative to a given timing solution. In 
each case we applied the given timing solution to the ob- 
servation times and added the phase residuals to obtain 
the phases as a function of time. We then fitted a tim - 
ing solution similar to the one used bv IStrohmaverl (|2004l ). 
i.e. (j){t) = 00 + i^{t - to) + u{t - to)^/2. We included an 
extra term ((jjp) because we fixed to to be the same as 
IStrohmaverl (|2004l '). so this value is no l onger arbitrary. We 
obtain the same fitted parameters as [Strohmayer (200^ 
but sl ightly different parameters to those of Ramsay et al.l 
(|2006t ). For reasons we sha ll explain, it is our fit to the 
data of lRamsav et al. I l|2006l ) that is given in Table \KT\ This 
corresponds to P = —3.21(61) x lO-^^^ss"^ which can be 
compared with the v alue P = —3.31(09) x lO-^^^ss-^^ from 



iRamsav et al.l l|2006l). Our unce rtainty is six times larger 
than that of Ramsav et al.l (|2006l ) and the values are slightly 
different because we allowed more freedom in the fit. We 
think our fit is the correct one because there is no reason 
that the fit has to have ze ro value and gradient at t = 0, as 
was effectively assumed bv IRamsav et all l|2006l ). who fitted 
only a parabolic term. 

Comparing the two ephemerides we found that the 
value of the frequency derivative was not consistent. More- 
over if we calculate the phase of the maximum of our ob- 
servation with the ephemeris from Strohmayer we obtain 
0.8170±0.0016 for the May 2003 and 0.7328± 0.0013 for Au- 
gust 2005 (different from 0.97 of Tablel6lfrom [Ramsav et al.l 
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I2OO6I ). This suggests that the two ephemerides have differ- 
ent zero points contrary to what was stated in the respective 
papers. Next we compared the phases predicted by the two 
ephemerides expecting to see a constant offset between the 
two. Instead we found a drift between one ephemeris and 
the other. The phase difference started at approximately 
zero for the first observation with ROSAT (so indeed the 
two ephemerides had exactly the same zero point) but were 
0.15 cycles apart for the last Chandra observation. Therefore 
we re-analysed the R OSAT data fr on i the 30"" April 1996 
and phase- folded it on IStro hmayer ''s (2004') ephemeris, but 
obtained the same phasing as Ramsav ct al. (2006). We can 
only obtain the same phasing as lStrohmaveil p004ll if we do 
not apply the UTC to TT correction. In the ROSAT doc- 
umentation it says that the times are in UTClj. The error 



in the ROSAT times causes IStrohmaveij (|2004f ) to under- 
estimate the rate of spin-up in V407 Vul, a nd this is why 
his fre quency derivative is lower than that of I Ramsav et al.l 
l|2006t l. 

To be sure of the correct phasing between the X-ray 
and optical light curves of V407 Vul we used Table Ell s 
ephemeris and applied it to Chandra data taken on 19"^ of 
February 2003 and 24**" November 2003, which were taken 
before and after our May 2003 observation. We obtained 
the same relative phasing of the opti cal and X-rays at eac h 
epoch. When we use the ephemeris of I Ramsav et al.l (|2006h . 
the two Chandra X-ray light cu rves are alm ost perfectly 
aligned (Figure [5]), but if we use IStrohmaveit s ephemeris, 
there is a distinct shift between them. We ta ke this as fur- 
ther evidence of a problem with IStrohmaveJ 's ephemeris. 

To conclude, we used the epheme ris that resulted from 
refitting the data lRamsav et ahl (|2006i l to give the ephemeris 
listed in Table lAll In the top panel of Figure lAll we 
show the residuals of our fitted phases for V407 Vul af- 
ter removal of a constant frequency model with i/q = 
0.0017562482721063 Hz. We also show the fitted parabola 
minus the linear fit. 

For V407 Vul we performed an F-ratio test for the 
parabola versus the linear fit. The value of the parabola 
is 23.4 and of the linear fit is 115, we have 10 points so 
we obtain an _F-ratio=4.33 which is significant at the 95% 
confidence level but not at 99%. 



AO. 2 HM Cue's ephemeris 

As was explained above in the case of V407 Vul, in order 
calculate the uncertainties in phase due to the ephemeris we 
need to know the covariance terms of the ephemeris. There- 
fore we applied the same metho d as before and dig itised and 
fitted the data of Figure 7 of IStrohmayer" ('2005'). We ob- 
tained the same fit coefficients as Strohmayor (2005) so our 
digitisation does not cause loss of information; in this case 
no ROSAT data were involved. We also obtain the covari- 
ance terms. Our fitted parameters are given in Table IA2I 
We show the phase residuals after subtracting a constant 
frequency model {vo = 0.0031101279743869 Hz) in the bot- 
tom panel of Figure lAll We applied an F-ratio test to the 
HM Cnc data. This time there were 69 points and we ob- 
tained a of 54.9 for the parabolic fit and 10380 for the 
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Figure Al. Phase residuals of V407 Vul (top panel) and HM Cnc 
(bottom panel) after subtraction of constant frequency models. 
The dashed line shows our parabola ephemeris subtracted the 
linear fit. 



to (TDB) 
<t>o 

i^(Hz) 
i>(Hzs-i) 
r(4>o,i') 
r{<t>o,i') 
r(f, i>) 



53009.889943753 

0.0003(14) 

0.00311013824(10) 

3.63(0.04) X 10-1^ 

-0.48041115 

-0.61096603 

0.94898169 



Table A2. Ephem eris of HM Cnc derived from Figure 7 of 
IStrohmaveJ ||2005| ) . The uncertainties of the parameter's are 
given within brackets. We also give the correlation coefficients 
for the fitted parameters. 



linear fit. This gives an _F-ratio of 186, significant at the 
99.99% confidence level. 



AO. 3 Uneertainties on phases from the ephemerides 

To calculate the uncertainties in the absolute phases due the 
uncertainties in the ephemeris we used the relation: 



L + - 



tQfo-l + {t~tofal/A + 



2(t ~ to)C^„^ + (i - to fCc + {t~ toYCc^„ (Al) 

where Cxy is the covariance of X and Y and can be written 
CxY = a-xo'Yr{X,Y). We give the correlation coefficients 
r{X,Y) in Table ED and Table IMl 

For the phase difference between two epochs Acp = 
0(^2) — 0(ii) one cannot simply combine in quadrature the 
uncertainties on the absolute phases at each epoch because 
the same coefficients are used in each case. (This is most 
easily seen by considering the case of two identical epochs 
for which the uncertainty in the phase difference must be 
zero.) Instead one must use the following relation: 
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We used this to calculate the uncertainties on the phase 
differences in sections 13.21 and 13.41 



